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METHOD FOR DETERMINING THE PHYSICO-CHEMICAL PROPERTIES OF A 

THREE-DIMENSIONAL BODY 

FIELD OF THE INVENTION 

The invention relates to determining the physico-chemical properties of a 
three-dimensional body; specifically, the invention relates to a method for 
determining the physico-chemical properties of a three-dimensional body. More 
specifically, the invention relates to a method for determining the mineral resources 
or reserves of a mineral body or layer. 

BACKGROUND OF THE INVENTION 

Several methods are known in the state of the art for determining the 
physico-chemical properties of three-dimensional bodies. Specifically, for 
determining the mineral resources or reserves of a mineral body or layer, this is, for 
calculating geological resources or mineral reserves in mineral bodies in the form of 
a layer. The most widely used methods are: 

Sections method: using bores made in sections that cut the mineral body, 
calculations are made obtaining the grades in each section. Then the area of each 
section is calculated and multiplied by half the distance to the anterior and posterior 
sections to thereby obtain the volume. Although the advantage of this method is that 
it can be applied to all types of layers, even very folded ones, it has many 
disadvantages, such as that each time a calculation parameter is changed, as the 
cut-off grade, the process must be started all over again; that as a grade calculation 
is made in each section, an interpolation direction cannot be used; that the bores 
not in the sections of calculation must be projected to the nearest one, complicating 
the process and, finally, that the sections method is very difficult to computerise. 

Polygons method: this method consists of projecting the centres of the 
intersections onto a plane and assigning to each intersection a polygon defined by 
the method of perpendicular or angular bisectors. Each polygon shall have the laws 
and powers or the intersection in the centre. Although this method is easy to apply 




and computerise, it has the following disadvantages: it cannot be used for folded 
layers; the calculation is not performed by interpolation of several bores, so that the 
grades obtained are over-optimistic; and it does not work in three dimensions. 

Triangles method: this method consists of projecting the intersections of the 
mineral layer onto a plane and defining the triangles formed when joining the 
vertices by triangulation. Each triangle is given the power and grades of the median 
of the intersections in the vertices. As with the previous method, this method is easy 
to use and computerise but it cannot be used for folded layers nor in three 
dimensions. 

Blocks method: this method consists of dividing the calculation area into 
blocks (parallelepipeds) and calculate the properties of each block interpolating with 
the intersections around it. This is the most widely used method, but its 
disadvantage is that for layer-shaped mineral bodies, as parallelepipeds are used, 
the geometric shape of the layer does not resemble the geometric shape of the 
blocks, and in thin layers it becomes even more complex. 

Thus, there is a need in the state of the art for an alternative method for 
determining the physico-chemical properties of a three-dimensional body that can 
improve on the commonly-used methods. 

The object of the present application is to provide an alternative method for 
determining the physico-chemical properties of a three-dimensional body, more 
specifically for determining the mineral resources or reserves of a mineral body or 
layer. 

The present method, which fulfils the requirements of working in three 
dimensions and being fully computerisable, is based on the iterative use of the 
triangulation method on the extrapolation of data obtained by bores. Moreover, the 
method of the invention illustrates the following advantages over the methods 
known in the state of the art: 

Any change of calculation parameter does not require a redefinition of the 

calculation units, 
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It defines calculation units in space, which can later be used to plan, draw 
and export to other programs, 

- It is possible to interpolate with any of the available methods, from the 
simplest method of assigning to each calculation unit the value of the 

5 nearest intersection, to applying the inverse of the distance or geostatistical 

methods. 

- It represents faithfully the power of the layer or mineral body, a 
fundamental information in thin layers. 

1 0 BRIEF DESCRIPTION OF THE FIGURES 



Figure 1 illustrates the drillings or bores made in a three-dimensional body or 

layer. 

Figure 2 illustrates the intersections resulting from the bore or drill segments 
15 that cut a mineral body or layer. 

Figure 3 shows a calculation unit, consisting of the part of the three- 
dimensional body or layer that has the same data (data 1, data 2) after the 
interpolation. 

Figure 4 illustrates the surface in space of the three-dimensional body or 
20 layer at its mid point defined by triangulation (Tl), this is, a set of triangles linked in 

space that define a surface in the centre of the three-dimensional body or layer. 

Figure 5 shows a cluster of points (NPS) generated by regular spacings in 
the two main directions of the three-dimensional body or layer. 

Figure 6 shows the new surface T2 (as well as a detail of this surface) 
25 defined by triangulation of the points of the cluster NPS. 

Figure 7 shows the three-dimensional representation obtained by applying 
the method of the present invention. 

Figure 8 illustrates the layer T1 defined by triangulation of the data obtained 
from the bores and their interpolation from example 1 . 
30 Figure 9 shows the cluster of points NPS and the surface T2 obtained by 

triangulation in example 1 . 

Figure 10 illustrates the golden grade of the mineral layer of example 1 . 



Finally, figure 1 1 shows the three-dimensional view of the mineral layer of 



4 



example 1 . 

DETAILED DESCRIPTION OF THE INVENTION 

5 To aid the comprehension of the present invention, the meaning of some of 

the concepts used in the present text is explained below: 

Three-dimensional body: a spatial body that may be predominantly in two of 
the three dimensions. When the method is applied to calculate geological 
10 resources, it will be a mineral body or layer. 

Bores: drillings made in bodies or layers to obtain samples for analysis and 
interpretation. 

Intersection: segment of the bore that cuts a layer of three-dimensional 

1 5 body. 

Interpolation: Calculation mode in which we define the data of a point of the 
layer or three-dimensional body using the information on the intersections 
surrounding it. It is possible to use the simplest method, in which the point is given 
the value of the nearest intersection, or the arithmetical mean of the intersections at 

20 a maximum distance, by an inverse power of the distance; or geostatistical 

interpolation methods, Kriging, etc.. It is also possible to use intersection search 
ellipsoids giving preferred directions, as is conventional in geostatistics. 

Calculation unit: part of the layer or three-dimensional body which for 
calculation purposes shall have the same Datal, Data2, etc. obtained from the 

25 interpolation. 

In a first aspect, the invention provides a method for determining the 
physico-chemical properties of a three-dimensional body that involves: 

30 a) Generating a database (BDS) that contains the data on the bores that define 

the situation and the physico-chemical properties of the three-dimensional 
body, 

b) Defining the surface (T1) in the spatial centre of the three-dimensional body by 
triangulation, 



c) Defining on T1 a cluster of points (NPS) generated with regular spacings in the 
two main directions of the three-dimensional body, 

d) Generating, by creating linked triangles between the points of NPS, a new 
surface (T2), very similar to T1 but in the suitable format for interpolation and 
graphical representation, 

e) Calculating, by any interpolation method, the properties of the points of NPS 
from the bore database BDS, 

f) Generating a new database (BDT2) using the triangles of the surface T2 that 
contains, for each triangle, the data of the coordinates of the vertices, the 
results of the interpolation of the vertices and the area of this triangle in space, 

g) Generating reports with the desired information from the database BDT2 and 

h) Generating three-dimensional graphical representations from the database 
BDT2. 

According to the method of the present invention, the database BDS is 
generated in stage a) from the information obtained in the intersections (see figures 
1 and 2) and comprises the following data: 

- Data n the (x, y, z) coordinates that define the position of each bore (s1 f 
s2, etc.) in the three-dimensional body (intersection of the bores and the 
three-dimensional body), where the coordinates can either define a 
single point that determines the centre of the body or an interval 
determining the beginning and the end of the three-dimensional body, 

- Data on the properties of the three-dimensional body such as the data on 
the actual width of the three-dimensional body (real power), analysis 
data, geotechnical data, geological data, etc. (data 1, data 2, etc.) for 
each bore (s1, s2, etc.). 

Then the stage b) is performed, in which the surface (T1) is generated in the 
spatial centre of the three-dimensional body by applying the triangulation method to 
the database BDS (see figure 4), specifically using: 

- The coordinates of the centre of the bores, 

- The three-dimensional interpretation of the known data of this body, 

- Prior knowledge of the typical shape of this type of body. 
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The triangulation method consists of forming linked triangles between the 
points that form the database. An algorithm is preferably used, such as the 
Delaunay algorithm. 

In the next stage, c), a cluster of points (NPS) is defined on the surface T1, 
generated by any algorithm based on regular spacings on the surface, this is, on the 
two main directions of the three-dimensional body (see figure 5). A possible 
algorithm can be as follows: 

- Generate the lines defining the intersection between the surface and 
equidistant parallel sections in each of the main planes, 

Divide these lines into equal segments, 

- The set of vertices defined by the lines in each segment shall form a 
cluster of points equidistant in one direction to the separation between 
the sections and in the other direction in the size of the segments. 

According to stage d), performing a triangulation on the points of the cluster 
of points NPS generates a new surface, T2, very similar to T1 but with the suitable 
format for interpolation and graphical representation (see figure 6). 

Then, in stage e) of the procedure the properties of the points of NPS are 
calculated by any interpolation method, ranging from the simplest method of giving 
it the properties of the nearest bore, a power of the inverse of the distance or any 
statistical method, using the bore database BDS. 

Then a new database is generated (BDT2) using the triangles of the 
previously generated surface T2 that contains, for each triangle, the data of the 
coordinates of the vertices, the results of the interpolation of the vertices, and the 
area of this triangle in space. 

Finally, the database BDT2 allows generating reports or graphical 
representations of the layer or three-dimensional body (see figure 7). Graphical 
software can be used to obtain the graphical representations, keeping in mind the 
following (see figure 3): 
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- Each triangle shall be the centre of a calculation unit, 

- Each triangle shall have in each vertex a segment that measures the real 
power at this point with the direction of the average of the perpendiculars 
to the planes formed by all the triangles sharing this vertex. In this way 

5 all triangles sharing a vertex also share this segment (edge) allowing all 

the calculation units to fit in perfectly in space, 

- The three aforementioned segments, together with the two triangles 
formed by joining their ends, define the volume of each calculation unit. 

10 A second aspect of the invention consists of applying the previously 

described method to determine the resources or mineral reserves of a mineral body 
or layer. This method comprises the following stages: 

a) Generating a database (BDS) that contains the data on the intersections of the 
15 bores defining the mineral body or layer, this database comprising: 

- Data of the (x,y,z) coordinates defining the position of each bore (s1, s2, 
etc.) in the mineral body or layer (the intersection of the bores with the 
mineral body or layer), wherein the coordinates can either define a single 
point determining the centre of the body or an interval determining the 

20 beginning and the end of the three-dimensional body, 

- Data on the properties of the mineral body or layer (data 1, data 2, etc.) 
for each bore (s1, s2, etc.). 

b) Defining the surface in the spatial centre of the mineral body or layer (T1) by 
forming linked triangles between the median points of each bore position (s1, 

25 s2, etc.) or intersections; to do so the following steps shall be followed: 

- Using the centres of the intersections of the bores with the mineral layer, 
the information on any outcrops of the layer and the geological 
interpretation regarding the spatial location of the layer, a set of points 
and lines are defined located on the central surface of the mineral body 

30 or layer, 

- Using these points and lines, the surface they form is defined by 
triangulation, providing a set of linked triangles in the space, 
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- As many points and lines are added so that the surface generated by 
triangulation is a faithful representation of the centre of the mineral layer 
or body and it covers the entire area to be included in the study; 

c) Defining on T1 a cluster of points (NPS) generated with regular spacings in the 
two main directions of the three-dimensional body, for which the following steps 
are followed: 

- An algorithm is used to fill in the surface T1 with points that are more or 
less equidistant to one another, 

- The distance between the points is defined according to the calculation 
detail required so that its final three-dimensional representation agrees 
with the initial interpretation of the layer, 

- Depending on the algorithm used, the real distance between the points is 
not necessarily always the same; 

d) Generating, by forming linked triangles between the points NPS, a new surface 
(T2) that will be very similar to T1 but has the suitable format for interpolation 
and graphical representation, for which a triangulation algorithm shall be used 
on this cluster of points, 

e) Calculating, by any interpolation method, the properties of the points NPS from 
the bore database BDS, 

- When interpolating, for each point of NPS the properties of the three- 
dimensional body at this point are calculated using the information on the 
intersections of the surrounding bores, 

- The interpolation can be by the simplest method of giving it the 
properties of the nearest intersection, a power of the inverse of the 
distance, or geostatistical methods such as Kriging or others, 

f) Generating a new database (BDT2), from the triangles of the surface T2, which 
contains, for each triangle, the data of the coordinates of the vertices, the results 
of the interpolation of the vertices and the area of this triangle in space, 

g) Generating reports with the desired information using the database BDT2. 

h) Generating a three-dimensional graphical representation from the database 
BDT2 by graphics software that allows a three-dimensional representation. 



In the same manner as described for the general method, when generating 
the three-dimensional graphical representation from the database BDT2 the 



9 



following shall be kept in mind: 



- Each triangle shall be the centre of a calculation unit, 

- Each triangle shall have in each vertex a segment that measures the real 
5 power at this point with the direction of the average of the perpendiculars 

to the planes formed by all the triangles sharing this vertex. In this way 
all triangles sharing a vertex also share this segment (edge) allowing all 
the calculation units to fit in perfectly in space, 

- The three aforementioned segments, together with the two triangies 
10 formed by joining their ends, define the volume of each calculation unit. 



The following example is allows illustrating the invention. 



EXAMPLE 1 

15 

A calculation is performed of gold (Au), silver (Ag), copper (Cu) and Arsenic 
(As) reserves of a mineral layer, specifically of the gold grade of this mineral layer. 
To do so, the following database is generated (BDS; table 1) from the data of the 
intersections of the mineral layer bores whose reserves are being calculated. 
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where: 

- (x1 f y1,z1) and (x2,y2,z2) are the initial and final coordinates of the 
intersection of the bore with the layer. 

- P_R is the real power of the layer in each intersection. 

- <Au>, <Ag>, <Cu> and <As> are the properties of the layer in each 
intersection, in this case they are analytical data of the elements Au, Ag, 
Cu and As. 

Based on the coordinates of the centres of the intersections and the 
geological interpretation, a surface (T1) is defined by triangulation that represents 
the centre of the layer (see figure 8). 

Then, the cluster of points (NPS) is defined on the anterior surface T1 
followed by the triangulation T2 (see figure 9). 

In this way, for each vertex we have its coordinates and the results of the 
interpolation, and for each triangle of T2 we have the information on the three 
vertices that define it, so that the triangle represented in the following table will be 
that formed by the vertices 30038000070, 30038500060 and 30039000060, where 
each vertex has real power (P_R) and <Au>, <Ag>, <Cu> and <As> values 
obtained from the interpolation of the intersections of the surrounding bores, which 
are also shown in the table. 

In this case the interpolation has been made by the inverse cube of the 
distance and the distances (Dist. in the table) are the distances between the point 
and the centres of the intersections of the bores. 

g = [ g, / ( d, ) p ]/[i/(d,) p ] 

g = result of the interpolation, 
g j = data of intersection i. 

dj = distance from the centre of intersection i and the point 
being interpolated. 
P = 3 
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The last row of the previous table represents the arithmetical mean of the 
P_R t <Au>, <Ag>, <Cu> and <As> values in the three vertices of this triangle, which 
together with the are of the triangle will complete all the information needed for this 
triangle when generating the reports with the calculations and for its three- 
dimensional graphical representation. 

Thus for example, separating in the database BDT2 the calculation units 
(triangles) with an <Au> grade over 4000 and grouping by categories, according to 
the nearest intersection, the following data table is obtained: 
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7504.29 


4.54 


4219.00 


12319.37 


4 


6940.00 


1.40 


7721.18 


6.02 


5625.34 


9479.52 


Total 


57017.00 


1.18 


8091.86 


7.74 


6638.76 


9605.54 



Figure 10 shows the triangles of the above tabie according to the <Au> 
grade. Finally, figure 1 1 shows a three-dimensional view of the calculation units 
generated with a 3D viewer. For the sake of a better three-dimensional 
representation the units have been slightly separated. 
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